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White dwarfs are almost completely degenerate objects that cannot obtain energy from 
thermonuclear sources, so their evolution is just a gravothermal cooling process. Recent 
improvements in the accuracy and precision of the luminosity function and in pulsational 
data of variable white dwarfs suggest that they are cooling faster than expected from 
conventional theory. In this contribution we show that the inclusion of an additional 
cooling term due to axions able to interact with electrons with a coupling constant <? ac ~ 
(2 — 7) x 10 -13 allows to fit better the observations. 

1 Introduction 

During the cooling process, white dwarfs experience some phases of pulsational instability 
powered by the k- and the convective-driven mechanisms pQ. Depending on the composition 
of the atmosphere, variable white dwarfs are known as DAV (atmospheres dominated by H) 
and DOV, DBV (atmospheres non dominated by H). These objects are experiencing g-modc 
non-radial pulsations, where the main restoring force is gravity. An important characteristic of 
these pulsations is that their period experiences a secular drift caused by the evolution of their 
temperature and radius. For a semi-qualitative purpose this drift can be well approximated 
by |5]: d\nU/dt ~ —adhiT/dt + bdlnR/dt, where a and b are constants of the order of unity 
that depend on the details of the model, and R and T are the stellar radius and the temperature 
at the region of period formation, respectively. This equation reflects the fact that, as the star 
cools down, the degeneracy of electrons increases, the buoyancy decreases and, as a consequence, 
the spectrum of pulsations gradually shifts to lower frequencies. At the same time, since the 
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star contracts, the radius decreases and the frequency tends to increase. In general, DAV and 
DBV stars are already so cool (and degenerate) that the radial term is negligible and the change 
of the period of pulsation can be directly related to the change in the core temperature of the 
star. Therefore, the measurement of such drifts provides an effective method to test the theory 
of white dwarf cooling. This, in turn, allows to obtain a simple relationship [3j [4] to estimate 
the influence of an additional sink of energy, axions for instance, on the period drift of variable 
white dwarfs: 



(ix/imodol) ~ (n o bs/n mo dcl) — 1 (1) 

where the suffix "model" refers to those models built using standard physics and Lx is the 
extra luminosity. 

Another way to test the theory of white dwarf cooling is based on their luminosity function. 
This function is defined as the number density of white dwarfs of a given luminosity per unit 
magnitude interval: 

rM u 

n{l)= / $(M)*{t)T cool (l,M)dM (2) 

J Mi 

where t satisfies the condition t = T — t coo i(l,M) — ipg(M) and I = — logtX/L©), M is the 
mass of the parent star (for convenience all white dwarfs are labeled with the mass of the 
main sequence progenitor), t coo \ is the cooling time down to luminosity I, t coo \ — dt/dM\, \ is 
the characteristic cooling time, M s and Mj are the maximum and the minimum masses of the 
main sequence stars able to produce a white dwarf of luminosity I, tpg is the lifetime of the 
progenitor of the white dwarf, and T is the age of the population under study. The remaining 
quantities, the initial mass function, 4>(M), and the star formation rate, ^(t), are not known a 
priori and depend on the properties of the stellar population under study. In order to compare 
theory with observations and since the total density of white dwarfs is not well known yet, 
the computed luminosity function is usually normalized to the bin with the smallest error bar, 
traditionally the one with I = 3. An important property of Eq. ^ is that the bright branch 
of the luminosity function is only sensitive to the average characteristic cooling time of white 
dwarfs at the corresponding luminosity when this function is normalized. 

pM u 

n = (r coo i) / <f){M)4>{T - i cool - t ps ) dM. (3) 

J Mi 

The reason [U [5] is that the stellar population is dominated by low-mass stars and, since the 
lifetime of stars increases very sharply when the mass decreases, the lower limit of the integral 
in Eq. ([3]) is almost independent of the luminosity, so the value of the integral is absorbed by 
the normalization constant. 



2 The case of G117— B15A and the luminosity function 

The measurement of the secular drift of the period of pulsation has been performed in the case 
of G117-B15A [B], a member of the ZZ Ceti (DAV) stars. The most recent value obtained so 
far is [TJ: 

(dIL/dt) ohs = (4.89 ± 0.53 ± 1.56) x 10" 15 s/s (4) 
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Figure 1: White dwarf luminosity function. The solid lines represent the models obtained with 
(up to down) g aee /10~ 13 = 0,2.2,4.5 respectively. 



with an estimated proper motion correction II = — (7.0±0.2)xl0~ 16 s/s. Theoretical predictions 
[8] indicate that this star should experience a secular drift of only II = 1.2 x 10 -15 s/s. Similar 
values (IT = 1.92 x 10~ 15 s/s or II — 2.98 x 10~ 15 s/s depending on the adopted mass of 
the envelope) have been independently obtained [9]. These results suggest that white dwarfs 
are cooling faster than expected (it important to confirm this statement by measuring this 
drift in other stars). There are three possible reasons for this, i) An observational error. 
This measurement is difficult and it has been obtained by only one team on just one star. 
Although the measurement tends to stabilize, it suffered strong fluctuations in the past [10] . 
ii) A modelling error. Models have been noticebly improved during the last ten years and the 
two independent models computed up to now [8j [9| are in a qualitative agreement, iii) An 
additional sink of energy is responsible of the accelerated cooling rate [3] . 

Under the conditions of temperature and density in the interior of G117-B15A, the axion 
emission rate is dominated by electron brcmsstrahlung (only the DFSZ axion model) that 
behaves as e ax oc ffi cc Tj erg/g/s, where Tj is the temperature in units of 10 7 K and g aeo is the 
strength of the axion-electron Yukawa coupling jllj . Thus it is possible to include the axion 
emissivity in Eq. (QJ and adjust g acc to fit the observed values. The value that best fits the 
observations is in the range of g acc ~ (3 — 7) x 10~ 13 . 

Figure [T] displays the white dwarf luminosity function, DA and non-DAs |12j . The values 
from the DR4 of the SDSS [13] (squares) were complemented at high luminosities with those 
obtained using the more recent DR7 [14] (crosses) . It is important to notice that models without 
axions [15] predict an excess of white dwarfs in the region \og(L/L Q ) ~ —2 as in the case of the 
pure-DA sample [3]. If axions are included, the best fit is obtained for g aco ~ (2 — 3) x 10~ 13 . 
The luminosity function at high luminosities is still poorly known from both the theoretical 
and observational point of views, but it is clear that will provide strong constrains. 
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3 Conclusions 



There arc two independent pieces of evidence, the luminosity function and the secular drift 
of DAV white dwarfs, that white dwarfs are cooling down more rapidly than expected. The 
introduction of an additional sink of energy linked to the interaction of electrons with a light 
boson (axion, ALP,. . . ) with an strength g aeo ~ (2 — 7) x 10~ 13 solves the problem satisfactorily. 
Naturally, the uncertainties that still remain, both observational and theoretical, still prevent 
to claim the existence of such interaction. A systematic analysis aimed to discard any possible 
conventional solution is under way. 
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